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Abstract 

Recently Eling and Oz [1] proposed a simple formula for the bulk viscosity of holo- 
graphic plasma. They argued that the formula is valid in the high temperature (near- 
conformal) regime, but is expected to break down at low temperatures. We point out 
that the formula is in perfect agreement with the previous computations of the bulk 
viscosity of the cascading plasma [2,3], as well as with the previous computations of 
the bulk viscosity of M = 2* plasma [4,5]. In the latter case it correctly reproduces 
the critical behaviour of the bulk viscosity in the vicinity of the critical point with the 
vanishing speed of sound. 

March 2011 



Contents 



1 Introduction and Summary 2 

2 Bulk viscosity of the cascading plasma 4 

2.1 Comparison of (1.2) with bulk viscosity of the cascading plasma .... 7 

3 Bulk viscosity of Af = 2* plasma 7 

3.1 Comparison of (1.2) with Af = 2* bulk viscosity away from criticality: 
5<5 C 9 

3.2 Comparison of (1.2) with Af = 2* bulk viscosity at criticality \5 C — 5\ <C 1 11 

A Coefficient V 4 13 



1 Introduction and Summary 

In [1] Eling and Oz (EO) considered 1 effective five-dimensional gravitational description 
of strongly coupled gauge theory plasma frequently arising both in phenomenological 
and the formal examples of the holographic gauge theory/string theory correspondence 
[7,8] 

/ <r>- x — - 2 }2 (^) 2 - v{<f>i) j + s„ . (i.i) 



s R = 



16G 5 % 

In the absence of chemical potentials for the conserved U(l) charges they proposed a 
simple formula for the ratio of bulk-to-shear viscosities in holographic plasma dual to 

where T is the plasma temperature (equivalently the Hawking temperature of the black 
brane dual to the equilibrium thermal state of the plasma), 

2 rf(lnT) 

c * = ^y (L3) 

is the speed of sound waves in plasma, and (j>f are the values of the gravitational scalar 
fields (pi evaluated at the black brane horizon. The remarkable feature of the expression 
(1.2) is in the fact that it can be evaluated entirely from the 'horizon data'. The reason 



1 See [6] for related earlier work. 
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why this is unexpected is that the bulk viscosity is sensitive to the microscopic scales 
in the theory, which in holographic correspondence are encoded in non-normalizable 
components of the appropriate scalar fields near the boundary 2 . In fact, above was 
precisely the reason the authors of [1] suggested that the formula (1.2) validity should 
be restricted to the high temperature case — in this case the scalar fields are expected 
not to flow a lot from the boundary to the horizon, and thus their horizon values 
should remain sensitive to the microscopic parameters of the boundary gauge theory. 
EO formula (1.2) was successfully verified [1] for a wide set of gauge theory/string 
theory examples, where bulk viscosity was evaluated from the sound waves attenuation 
coefficient [14-17]. Furthermore, it was found [1] that while (1.2) is valid in some 
phenomenological models of gauge/gravity correspondence [18-20], it is violated in 
some other models [21,22]. 

In this paper we demonstrate that, at least in the context of formal examples of 
holographic correspondence, the range of validity of (1.2) is much wider. In particular, 
in section 2 we verify (1.2) for the cascading gauge theory plasma to order O (in -4 ^) 
in the high-temperature T ^> A expansion 3 . Results are presented in section 2.1. In 
section 3 we verify (1.2) for M = 2* gauge theory plasma for all temperatures 4 . In 
particular, we demonstrate that (1.2) is valid in the vicinity of the phase transition 
with vanishing speed of sound [23-25]. Notice that in this latter case the finiteness 
of the bulk viscosity at criticality, along with the vanishing of the speed of sound at 
criticality, implies that the derivatives of the scalar fields must diverge at the critical 
point for (1.2) to be valid. Results are presented in sections 3.1 and 3.2. 

It would be interesting to understand precisely why (1.2) fails for phenomenological 
models [20-22] . It is intriguing that these models also violate the bulk viscosity bound 
proposed in [4]. Thus, the derivation of (1.2) might be the first step towards proving 
the bulk viscosity bound in formal holographic examples. 

2 This should be contrasted with the familiar universality of the shear viscosity to the entropy 
density ratio [9-12] which can be evaluated in the membrane paradigm framework [13]. 
3 Original computations of the bulk viscosity in this theory were done in [2,3]. 
4 Original computations of the bulk viscosity in this theory were done in [4,5]. 
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2 Bulk viscosity of the cascading plasma 



Effective gravitational action for the cascading gauge theory plasma was derived in [2] 

* ■ as /» {»-"(«/>'- mm' - i(»>' - d." >'•--"" 

-'}• 

(2.1) 

where 

7> = -24e~^ f ~ 2w + 4e~~ f ~ 12w + p 2 e ^'~ f+4w + -ATVt' . (2.2) 

2 

In this case (1.2) reads 5 



EO 



d$ H \ 2 80 (df H V (dw H \\ 1 ~ H _ Ar _ iwH (dK H ^ 2 



(2.3) 



As argued in [3], it is technically challenging to study bulk viscosity of the cascading 
plasma for all temperatures. Instead, it was studied in [3] perturbatively (to the third 
sub-leading order) at high temperature. From (2.3)-(2.7) of [3] 

f 2 {x) = d + S J2 { (^") (- 2 6n(^) + Wz) + ^A 2n (x)^ | , (2.5) 
f 3 (x) = a + 5 ^ | (-2£ 2n (:r) + W*) - ^A 2n (x)^ | , (2.6) 

= + E{ «*,(*)} , (2.7) 



where 



/= -ln/i + -ln/ 3 + — ln/ 2 , w = — In fy , $ = In g \ . (2.9) 
4 5 J 10 J ' 10/ 2 f 



3 We use so to distinguish the expressions obtained applying (1.2). 
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Notice that dependence disappears for the scalars entering (2.3), and the only tem- 
perature dependence enters via K± = K+(T). We refer the reader to [3] for the asymp- 
totic parametrization of various functions {/t 2 „, £,2n, V2n, A 2n , (2n}- 

The pressure V, the energy density S, and the entropy density s of the cascading 
plasma are given by [26] 



V 3 f 7 \ £ 3 / 4 \ /o . 



with (see (2.58) of [3]) 



7 P 2 /7 20 35 7 >20 7 , \ P 4 / 35 , n 35 20 7 , „ 
a 20 = + -k 2 ' + — C? + — n2 J — + hi 2 nl fi + — n 2 2 

+ 1Z!> + 1^,0 ln 2 + l^,o j 2 _7_ C 2 >0 7^,o _ 35 20 \ [ 175 

108 24 S1 12 2 12 S2 6 3 36 S1 ) \72 

7 1.9 n 175 9 n 35 9 n 1 ~ 35 , 9 n 1 ~ 875 175 ,9 n 35 _ 9 „ 

+ —Q + «V «V In 2 Ci In 2 + Cf In 2 

12 S3 54 2 18 2 36 S1 324 108 Sl 48 

+ L 1„3 2 + _^J,o ln 2 2 + _L^o ln 2 + 7_ 4 o ln 2 2 + ^2,0 ln 2 + 7 K j,o 

36 C2 18 3 )Ky U \K*. 

(2.11) 

The precise temperature dependence of was determined in [27] 

1, /64vr 4 sT\ 
_i = _ln( x — . (2.12) 

P 2 2 V 81 A V 

Using (2.12) and the expressions for the pressure and the energy density from (2.10) 
we find 

c2 = &P = 1 7 - 12a 2 , - 6P 2 §ff 
s ^ 3 7 + 4a 2i o + 2P^ ' 

Thus, given the perturbative high temperature expansion for a 2j0 we can evaluate from 
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(2.13) the perturbative high temperature expansion for c 



2 4P 2 /10 2, n 8 20 4 /20 \ P 4 / 2 a20i n 

3C - = = 1 ~ 3 IT. + ( Y " 3 2 " 3 K " - 3 C > ) 5f + (" 8 " 3 C? 2 

+ — In 2 - i In 2 2 - -Kj fl + _ 1 4° ln 2 + ^ _ 4 2 ,„\ 

3 3 3 3 9 Sl 3 2 9 2 3 S2 / AT| 

+ (4ln 2 2C?'° - ^ 2 - ^4°ln 2 2 - |cJ'°ln2 + ^(4°) 2 + y«J°C?'° 
+ l( Cl 2 ' ) 2 - i In 3 2 - 12 In 2 + | C 2 '° In 2 + ^ + ^4° In 2 + | In 2 2 



212 2 ,o 106 2>0 46 2j0 92 2>0 4 2>0 8 2>0 \P 



n K 2 n S>1 t- C, 2 + n K 3 0^3 Q K 4 7^4 + U 



* 



>io 



(2.14) 

Consistency of the first law of thermodynamics (which was verified in [3] with an 
accuracy of oc 10 -7 ) implies (see (4.5) of 

dK+ 



dlnT 



2 P 2 + [ ~ + 2d 2,0 + 4«|° + ln 2 



M + V 3 3 2 5 (Al ' J 3 Cl + 18 



60(77?, h ) 2 + 2C 2 2 '° + 2kY ln 2 + 1 In 2 2 - 3 ln 2 + ln 2 C 2 '° - 24£ 2 ^ 



2,0 



+ (f In 2 + g C ?'° + f 4° - \CI' Q - y 4° + 2Cf + 4. 2 '° + 16& - I In 2 Cl 2 

- 74° ln 2 - ^ In 2 2 + ln 2 C 2 2 '° + - In 2 2 Ci'° + 2 ln 2 4° + In 2 2 4° + - In 3 2 

+ 40«,) 2 + ^(A%) 2 - 36£% - 180*% vl h + y»h% (A%) 2 - yA% A% 

24 4 

- 240(*%) 3 + -(A%) 3 + 180^% V% ~ 30 In 2 (*%) 2 - 12 In 2 £ 2 % - - In 2 (A%) 2 

- 607 - 120(,%) 2 4° - 60(,%) 2 C 2 '° - 48^% 4° - 24^,, 4° - ^(A%) 2 4° 



108 



(2.15) 
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We now have all the necessary ingredients of evaluate (2.3) to order O {i^J 

c 

v 



EO 



8 P 2 ( 76 16 >20 32 20 8, n 8 >n 8 n 16 , n \ P 4 
= 9 K. + " 27 + ¥ C ' + T 4 + 9 2 - 9 C? '» + 3 + 45 A ° * « 



,.2,0 
2 



+ (l2 + 1 4° In 2 + ^ In 2 Cl 2 '° + f 4° C 2 '° - jjtf, A - ^C?, A - f Cf> 
+ y < fc Cf° + f + | A?,,C 2 '° + | A? >fc 4° " jjtf,* ^ 2 + f < h In 2 

+ ^ln2 - ^tf >fcC ?'° + |<A " 8 Ci'° " + f C?, fc - f < h - 

g"^2,fe 2"^°,^ + 45 A°,fe — g^2,/i + ~g~^ 2 ' ~g K ' 3 ' 3^2,71 27" ln2 + -(Ci^) 2 
+ jj(ln2) 2 + I (4°) 2 + jj(tf ) 2 - l f{vl h ? - |^(A?,J 2 ) ^ 

+ v 4 — + o 



(2.16) 

where the coefficient V4 is given in Appendix. 

Eq. (2.16) should be compared with the expression for the bulk viscosity derived 
in [3] (see (4.14) in 



C 8 P 2 4 / P 4 P 6 P 8 \ ( P 10 \ 

h^ + ! r 2 + fe + fe4 *?) + ° { -m ) ■ (2AT) 

For convenience we collect in Table 1 numerical values for all the coefficients entering 
(2.16) and (2.17) as computed in [3]. 

2.1 Comparison of (1.2) with bulk viscosity of the cascading plasma 

Using (2.16) and (2.17) and the data from Table 1 we find: 

C 8 P 2 P 4 P 6 P 8 / P 10 

- = + .1763445167— - 2.263612386— + 3.026509967— + O — 

V eo 9^ K 2 K% Kt \K* 

C 8 P 2 P 4 P 6 P 8 /P 10 

- =-— + .1763444983— - 2.2636127883— + 3.026511143— + 01 — 
V 9fT* Kl Kl K$ \Kl 

(2.18) 

3 Bulk viscosity of J\f = 2* plasma 

Effective gravitational action for Af = 2* gauge theory plasma was obtained in [28] 

1 



4ttG 



/ ^ 5 v^ [\R - 3(&0 2 - (9 X ) 2 - V] , (3.1) 

5 J Mr, 



n 


1 
1 


_ 


Q 
O 


A 


9 




0.73675974 


-0.62226255 


-0.03784377 


d n 

In 


-0. 01717287 


0.00534036 


-0.01064222 




A ' 


-0.87235794 


-1.11562943 


1.39008636 




Sn 


-0.15342641 


0.62226267 


-0.32514260 




o 

K n,h 




0.62226259 


-0.42061461 


0.00816831 


cQ 




-0.07981931 


0.01661150 


-0.00920379 


£0 




0.01919989 


-0.05277626 


0.01385333 


o 


-0.14891337 


-0.21809464 


0.00213345 




Ki,h 


0.16806881 


-0.14619173 


0.01639579 




AO 
Sn,/i 


-0.41123352 


0.33024116 


-0.07445122 




fan 




0.13225837 


-1.69770959 


2.26988336 



Table 1: Coefficients entering (2.16) and (2.17). 
where the potential 6 



V = — 



1 

16 



i r aw 

3 V da 



+ 



dW 

dx 



-W 2 

3 



(3.2) 



is a function of a and \i an d is determined by the superpotential 

1 



W = 



-2a 



e 4a cosh(2x) • 



In this case (1.2) reads: 

c 

v 



-ctT 2 



EO 



24 



da 



H 



dT 



+ 



H 



dT 



(3.3) 



(3.4) 



In what follows we present results only for Af = 2* gauge theory plasma at vanishing 
fermionic mass 7 [23]. The latter case corresponds to setting x = 0. 

As pointed out in [23], and further explored in [24,25], strongly coupled Af = 2* 
gauge theory plasma undergoes a second-order phase transition at 

5= (^±y = 5 c = 5.4098(6). (3.5) 



6 We set the five-dimensional gauged supergravity coupling to one. This corresponds to setting the 
radius L of the five-dimensional sphere in the undcformcd metric to 2. 

7 Although we did not verify this fact explicitly, based on the agreement for m/ = 0, we do not 
expect discrepancy between (3.4) and the direct computation of the bulk viscosity for m/ ^ 0. 



8 



1.00 



0.95 



0.90 



0.85 




Figure 1: (Colour online) Horizon value of the scalar field po = ^ as the function of 
5 = (^fr) 2 - Data points are blue; the solid red line is the best fit of the data with order 
15 polynomial in 5. The dashed vertical green line represents the critical point of the 
theory 5 = 5 C . 



At the critical point the speed of sound vanishes. We discuss separately the cases 
5 < 5 C , and \S C - 5\ < 1. 

We use asymptotic parametrization of the scalar function, the speed of sound, and 
the bulk viscosity as in [5] 

ft C 4 



p = lna, c s = — L, - = ^(/3 2 -l). 

Vo Tj 6 



(3.6) 



3.1 Comparison of (1.2) with J\f = 2* bulk viscosity away from criticality: 

5 < 5 C 

We rewrite (3.4) assuming functional dependence p H = liaa 11 = po(5) 



326 2 Pt fdpo 



From eq. (2.11) of [5], 



EO 3 P0 \ dS 



24vr 2 6 
o = — — pne D 



y/2 



(3.7) 



(3. 



Thus, given the data sets {pn, pio, Po ; a o, a i, Pi, P2} obtained in [4,5] we can reconstruct 
the functional dependence Pq(5). The result of such reconstruction is presented in 
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Figure 2: (Colour online) Comparison of the EO prediction for M = 2* plasma bulk 
viscosity with the explicit computations from the quasinormal modes [4]. The dashed 
vertical green line represents the critical point of the theory 5 = 8 C . 

Figure 1. Data points are blue — there are altogether 320 points. The solid red line is 
the best fit to the data with the polynomial 



15 



fit 

Po 



(3.9) 



i=0 



Having an (approximate) analytic expression for po ~ Pa ($), we can compute the 
prediction for the bulk viscosity of M = 2* plasma from (1.2) 



EO 



325 2 f3f ( dp, 

^pT 



fit 



d,5 



(3.10) 



We find that the numerical agreement between (3.10) and the original result (3.6) is 
excellent — Figure 2 presents the functional dependence of 



CM 



EO 



C/v 



l 



5. 



(3.11) 



Again, the vertical green dashed line represents the critical point of M = 2* plasma. 
The agreement is relatively worse for small values of S and in the vicinity of the critical 
point. In the former case, this is caused by numerical errors for the evaluation of the 
bulk viscosity ( the larger errors are induced due to small amplitude profiles of the 
gravitational scalar field a ), and in the latter case the relatively worse agreement is 
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0.83325 
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Figure 3: (Colour online) Horizon value of the scalar field p = e ° an d 5 = \pr) as 
a function of ft = c 2 at criticality. Data points are blue; the solid red lines are the best 
fit of the data with polynomials in ft — see (3.14). 

caused by the choice of the 'fit' (3.9) (in the vicinity of the critical point ^ is expected 
to diverge as /3f 2 ). As was already pointed out in [1], the agreement between (3.7) 
and (3.6) can be established analytically to order 0(5). In the next subsection, with 
a suitable parametrization of the horizon value of the scalar field at criticality, we 
drastically improve the agreement between (3.7) and (3.6) for 5 close to 5 C . 



3.2 Comparison of (1.2) with J\f = 2* bulk viscosity at criticality \5 C — 5\ <C 1 

An important prediction of [4, 5] is that the bulk viscosity of A^ = 2* plasma remains 
finite at criticality, where c 2 = 0. The EO expression (3.7), if correct, implies that 
the derivative of the scalar field po with respect to 6 must diverge as c~ 2 at criticality. 
Thus 5 parametrization of po is n °t very useful. Instead, in the vicinity of the critical 
point we rewrite (3.7) function of 



ft^cl 



ftl 
3 



We find 



EO 



965 2 



dpp 

dft 



1 d5 



(3.12) 



(3.13) 



where now we understand po = Po(ft) an d 5 = S(ft). These two functions can easily 
be reconstructed from the data sets {pn, pio, Po, do, ai, fti, fti} obtained in [4,5] — see 
Figure 3. Blue dots represent the data points, the solid red lines are the fits to the 
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Figure 4: (Colour online) Comparison of the EO prediction for Af = 2* plasma bulk vis- 
cosity with the explicit computations from the quasinormal modes [4,5] in the vicinity 
of the critical point, i.e., , for (3 = 0. 



data: 

5 fit =5.40987 - 1.14476 (5 2 - 4.64346 /3 3 - 14.55312 /3 4 , 
p f a =0.83305 + 0.09281 + 0.17683 (3 2 + 0.36473 /3 3 . 

We can now compute the prediction for the bulk viscosity of M = 2* plasma from (1.2) 



(3.14) 



96S 2 / dp J 



fit 



Pi 



eo /'(i \ dp 
Figure 4 presents the functional dependence of 

XM\eo 







d(3 

dW 1 



- 1 



vs 



(3.15) 



(3.16) 



in the vicinity of the critical point of M = 2* plasma. 
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A Coefficient V4 

v _ 308 V,o , 6160 2,0 _ 404 244 808 o 64 608 _ 64 

|/ 4- gl Si + gl ^2 27 ^+ g Vl,h+ 135 \h+ g K 2,h + g $2,h ^ ^ 

32 _ 80 2 ,o _ 160 ,0 _ 16 o _ 32 8 16 _8 

32 2,0 , 16 2,0 32 64 2j o 128 2 ,o 16 2 ,o 

+ ~g~ K 4 + ~g~~S3 ^Sl,/i 7 7l,h A l,h + ^Vl,h A l,hSl + ~^Vl,h A l,h K 2 ~^Sl,hVl,h^l 

32 ^° n° ,<- 2 ' 32 ^° X° ^ 2 '° 64 ^° X° ,<- 2 ' 32 ^° 2,0.2,0 , 32 ,.2,0.2,0 
~^^l,h'll,h K 2 ^gSl,/i A l,hSl 4gSl,h A l,/i re 2 ~g~Sl,h K 2 SI "T ~^'ll,h K 2 Si 

6_4 o 2,0.2,0 _ 272 2,o>2,o 32 64 o 512 2 ,o _ 176 ^2,0 

45 ,h 9 3 ^ h ' 1 ' h ^^h A i,h ^ 27 ^ 1 >' 1 2 9 'i.^ 1 

352 «° „-2,o _ 512 a2; o 1024 ^ 2 256 2)0 _ 128 

g '/i,/^ 135 135 i.h 2 27 45 'hh A i,h 

^ Vll,h) A l,h^ j^'ll,h\ A l,h) + g Sl,/iS2,h g l'/l,/J ^2 g l'/l,/J Si 

256 2,0 128 .2,0 H84 m0 2 2 ,o 592 2 20 112 2 
"g - S2,/i^2 _ ^ _ S2,hsi 225 ' 225 ' ~^ g - si,/A'/i,/J 

+ (\° ^ 2 - -t° n° - —C° A - -n° r° 4- -n° (7° ^ 2 - —A C° 

225 ' ' 3 < »l,ft''2,/i 4gSl,/i^2,/i 2'/l,hS2,^ ^ g'/l,fcVSl,/J ^ A l,AS2,/i 

+ 7° ^ 2 - — K° f 2 >° - — K° K 2 '° + ^/^O 2,0 4- I^V 2 ,0.2,0 4- ^1k 2 '°K 2 '° 

45 iMsi./J g "-2,/iSl g K 2,h K 2 + g Si ^3 + g Si S2 ^g^^S 

32 2,0 .2,0 _ 16^o ^2,o_ 32 2,0 16,. n2.2,o 32 \2 2,0 , i6 ^2,0 

~r "g" 2 S2 ~ ~g~~S2,ftSl "g~S2,/l K 2 ~g~VSl,jJ Si T "g~VSl,/J K 2 "T "^"'/2,/iSl 

32 20 32 >2 ^4 . 20 32 20 l^.o .20 ^*4^o .0 

+ ~~^~V2,h K '2 + ^ A 2,hQl + ~^ A 2,h K 2 g"Cl,/i /c 3' g _ Cl,/iC2' + ~g~Sl,/iS2,fc 

4. k° - (T 2 ' ^ 2 - — r° fK 2,0 ^ 2 4- —7i° k 2 ' 4- — n° C 2fi - —n° k° 

+ g ^l,h K 2,h gSl,?iVSl ) g Sl,/A K 2 J + g '/l,/i K 3 + g '/l,hS2 g r ll,h K 2,h 

+ ^?,,(Ci 2 '°) 2 + f <,(4°) 2 + |A?^ 2 ' + fA?^ 2 ' - f A?, A 

32 304 3296 28 2 272 2>0 68 2 , 2 1364 

45 1 ' h " 2 . ft ~y^~ A l,h A 2,h g "g"VSl,W "g - 1 K 2 J ej^ 1 ' ~^ g - V r /l,?iJ 

1064 2672 352 8 ,3 , / 436 2>0 872 2>0 

"225~ iAl ' J -9—^1,^ + n25 iAl ^ 9^1.^ + ^ "27~ K2 

-Ci\ - 16,?, - *?A? fc - — ^ h + — A°ft - -C 2 °, + ^Cf + ^4° 
27 135 9 3 1 45 ' 9 ' 9 9 

8 39 - - 8 1fi 39 Ifi 39 



+ 
+ 
+ 
+ 
+ 



27 135 L ' n g z > ra 3 ^^>' t 45 z > re g ,z > ra g ,z 9 

8 „o 4. 3_2_ 2,0.2,0 _ 8 _ ^6 A o _ 3j^ 2,0 , 16 ^2,0 32 K 2,o 

3 'l2,h^ g ^2 Si ^Sl,h'll,h ^Sl,h A l,h g Sl,,^ ^ 3 '/ljftSi ^ 3 '11^2 

32 .2,0 64 2,0 16^o .2,0 32 , 2378 8 2 32 , 2,0^2 

^ A l,fcSl -!-^g^l,h«2 "g-Sl^Sl "I" ^%h A l,h + -gj~ + glSl,hj + Y^2 J 

-fC 2 'V-— fr7 l 2 -^ ) 2 Vn2+^C 2 ' - — + -^° + 2r ? ° + —A 
gisi J 3 225 l A 1; J 27 +3^2 +^i,h+ 15 Ai )fc 

9 \ /i 



-|c?.)ln 2 2 + ll„32 



(A.l) 
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